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We investigate the limits of intracavity absorption spectroscopy with nonlinear media. Using a 
common theoretical framework, we compare the detection of a trace gas within an undriven cavity 
with gain near and above threshold, a driven cavity with gain kept just below threshold, and a cavity 
driven close to the saturation point of a saturable absorber. These phase-transition-based metrology 
methods are typically quantum-limited by spontaneous emission, and we compare them to the empty 
cavity shot noise-limited case. Although the fundamental limits achievable with nonlinear media do 
not surpass the empty cavity limits, we show that nonlinear methods are more robust against certain 
technical noise models. This recognition may have applications in spectrometer design for devices 
operating in non-ideal field environments. 



I. INTRODUCTION 

Trace gas identification methods play an important role in precision optical metrology with multiple applications 
having potentially large societal impact. These include explosives detection [I], manufacturing process control, and 
medical diagnostics, such as detection of bio-markers in human breath [2]. Many comprehensive reviews have compared 
various trace detection techniques, considering both practical implementation pQ as well as fundamental quantum 
limits [3|. One established technique is absorption detection of a trace gas within an optical cavity. A laser beam 
interacts with the absorber multiple times via high-reflectivity cavity mirrors, which effectively enhances the path 
length and amplifies the minute single pass absorption signal for a fixed volume of gas. 

Practical implementations of cavity enhanced detection are aided by the technique of cavity ring-down spectroscopy 
(CRDS), whereby the absorption is extracted from a fit of the decaying cavity output [4 . CRDS is largely immune 
to input laser intensity noise, and its practical utility has led to its prevalence in a variety of commercial applications. 
Furthermore, this technique can be made highly parallel by feeding the cavity with multiple frequencies resonant 
with adjacent cavity modes from a femto-second laser source [2]. Modern mirror and cavity technology achieves 
finesses of approximately 10 5 for IR wavelengths with corresponding resolvable absorption coefficients in the range of 
10 -14 cm -1 in a one second timescale [3]. This sensitivity achieved under optimum laboratory conditions represents 
the comparative benchmark for novel approaches. 

However, when considering practical field- deployed devices, background technical noise may significantly degrade 
performance. Therefore, in comparing absorption sensitivity across different technology platforms, one must consider 
robustness to non-ideal environments. Measurement schemes with nonlinear media, especially when coupled with 
modulation, have shown great promise in other fields for providing robustness to such technical fluctuations. For 
example, transistors, Josephson bifurcation amplifiers [5], and fiuxgate magnetometers [6 , all exhibit robustness 
in field- deployment. It is thus natural to consider the implications of inserting nonlinear media into intracavity 
spectroscopy systems. 

The simplest addition is an amplifier which can compensate the cavity loss, thereby increasing the effective cavity 
finesse. This approach is limited by spontaneous emission as is described via amplifier uncertainty relations in Ref. [7]. 
This limit can also be cast into the language of the 'quantum no-cloning' [8] and the fluctuation-dissipation theorems 
[9|. The gain element does not need to be driven externally, and a laser cavity (seeded by spontaneous emission) can 
also be used for sensitive trace-gas absorption spectroscopy. This was recognized early in the development of laser 
physics [10] and work on broad-gain medium lasers with relatively narrow trace gas absorption lines continues [11] . In 
Ref. [12] , Kimble addresses the specific example of an intracavity, narrow gain- medium laser without external drive, 
motivated by the observation that a laser has a nonlinear threshold point which is sensitive to intracavity absorption. 
In the absence of spontaneous emission, the laser output power curve exhibits a first order phase-transition-like feature 
at threshold. However, as recognized in Ref. [12], spontaneous emission smooths out the sharpness of the threshold 
point and limits the sensitivity. Nevertheless, for absorption measurements in the presence of technical noise, enhanced 
performance is still possible. 
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This idea of enhancing the absorption sensitivity by using the trace gas as a perturbation which kicks the system 
past a nonlinear bifurcation point (or phase transition) can also be applied to situations without gain. For example, 
optical bistability can occur for a saturable absorber within a cavity driven near saturation p~3j[14]. The introduction 
of additional absorption kicks the system past saturation, thereby amplifying the signal. While optical bistability is 
typically discussed as a potential nonlinear transistor element for optical computing, it can also be used for absorption 
spectroscopy; yet, the performance is once again limited by spontaneous emission. 

In this paper we review absorption spectroscopy techniques using nonlinear elements from a unified perspective and 
derive both fundamental and technical limits. The potential benefit of such nonlinear techniques depends critically 
on technology limitations, such as mirror coatings and light sources as well as on technical noise of a non-fundamental 
nature. We use a common stochastic formalism to analyze different regions of the allowable phase space with drive and 
gain as the independent parameters. In each case, the derived limits are compared to the empty cavity benchmark. 
We express the absorption coefficient uncertainty as a function of time, recognizing that measurement noise is not 
purely white, and dynamical timescales are important. We conclude by demonstrating that nonlinear techniques can 
approach, but not improve upon, the fundamental sensitivity limits of empty-cavity methods. Most importantly, we 
show that nonlinear approaches can robustly outperform linear counterparts in the presence of technical noise, as 
might be found when operating outside of a laboratory environment. 

Although we will focus on sensitivity under fundamental and technical limits, one important issue we will not treat 
is selectivity, whereby signal shapes are interpreted in time and frequency space to determine what species or mix 
of species is present. There are several ways to imagine incorporating the systems discussed in this paper into a 
selective apparatus, for example using diffusion metrics with a cavity gradiometer or femtosecond frequency combs 
with a broadband gain medium. 

The structure of this paper is as follows. We begin by framing the problem and defining the estimation terminology. 
Next, cavity terminology is reviewed and the optical shotnoise-limited empty cavity benchmark is presented in these 
estimation terms. We then introduce the dynamical modeling which underlies the nonlinear treatment, using a 
stochastic formalism to represent spontaneous emission. Three different regions of phase space are then treated in 
turn: gain without drive, gain with drive, and saturable loss with drive. In the latter case, we focus on nonlinear 
modulation techniques and compare to CRDS. For each case, the parameter range for the nonlinear technique with 
increased robustness in the presence of background noise models is demonstrated. 



II. ESTIMATION THEORY 



In order to treat both the linear and nonlinear absorption techniques with a common approach that renders 
comparisons meaningful in the presence of finite bandwidth noise, we first define our estimation terminology. An 
absorption spectrometer produces a spectral plot where, at each wavelength of input probe light, the output is used 
to determine the absorption coefficient of the sample, denoted as, in units of cm -1 . The precision of a single shot 
measurement as a function of time (at a particular input wavelength) can be represented by the Allan variance 
[15]. For a scenario with additive white noise only, the Allan variance initially integrates down as 1/T where T is 
the integration time. The pre-factor on this portion is termed the 'short term sensitivity' and is given in units of 
cm~ 1 /\ /r Rz and is typically due to white noise sources. Repeating uncorrelated measurements in time and adding 
noise in quadrature also results in a noisefloor with the same units. Due to slower noise sources the Allan variance 
typically levels off at a level termed the 'detection sensitivity' (or bias stability) and is represented in units of cm -1 . 
The sensitivity can be inferred from the short term sensitivity and the stated averaging time. 

We estimate the quantity x through the measurement y with the two being related by a function y = f(x,pi) which 
is inverted to give x = g(y,Pi). Here pi are generic parameters that may be noisy, and thus add to the overall error. 
The propagated uncertainty estimation is usually given by adding the uncorrelated noises in quadrature as 

(1 , 

With y as the measurement and x as the desired variable, related by y = f(x,pi), it is typical to refer to the quantity 

as the responsivity of the detector. In photodetectors, e.g., this is given in units of Amps per Watt, A/W. The 
normalized responsivity is defined as 

- r (df/dx) 

R= J = — f— (3) 



3 



We subsequently do not use R and refer to R as the responsivity with the implicit understanding that it is normalized. 
The sensitivity of the detector dx 2 can be represented (ignoring pi terms above) as 

dg\ 2 _ dy 2 _ dy 2 _ dy 2 1 



Since most experiments display a dominant fast noise that must be integrated down until a technical noise limit is 
met, the typical form of the uncertainty as a function of time (i.e., the Allan variance) is 

where 7 is the fluctuation rate of the noise source (in frequency space the noise being white below 7 and rolling off 
above) and vo is the short time variance of the signal, which may depend on 7 in certain circumstances. In the following 
analysis, can either be due to shotnoise, with the bandwidth determined by the decay of the cavity and inversely 
related to the intensity, or it can be representative of near-threshold laser dynamics which have larger variances and 
slower dynamics (but smoothly merges into shotnoise far above threshold). The variance vt represents a technical 
noise source, such as alignment or gain instability in the post-cavity detection. This noise could be represented as 
vt/0- + 7t£)> but we consider this bandwidth 7^ as too slow to be integrated past; hence the denominator is taken 
as unity. 

Generally, non- modulated schemes as presented below take the form of Eq. ([5|. When one ignores vt, there is 
typically no enhancement in the fundamentally limited sensitivity achieved by adding an extra nonlinear element. 
However, the nonlinear element can enhance the responsivity R such that when vt is non-zero, the long term perfor- 
mance improves significantly. In other words, the nonlinear system is not optimal in a fundamental sense, but makes 
the system robust to certain types of technical noise. 

The above propagation-of-uncertainty style of analysis is typical but not generally optimal. Instead, we could apply 
an optimal Bayesian estimator to the problem, utilizing full knowledge of the system parameters [26 . In the simple 
case where the noise and distributions are gaussian and the dynamics linear (or appropriately linearized), the solution 
for the uncertainty as a function of time (i.e., the Riccati equation solution) typically takes the intuitive form of the 
estimation variance in Eq. Nevertheless, our analysis suffices to evaluate the relative merit between linear and 
nonlinear absorption techniques. 

III. SINGLE PASS AND EMPTY CAVITY BENCHMARKS 

The simplest absorption measurement, and therefore the most transparent application of the above analysis, uses 
input power Pq incident on a medium such that the output power Pi is given by 

Pi = exp((a G -a L - a s )L)P (6) 

where ac is the gain per unit length, is the non-sample loss, as is the sample loss we are trying to detect, and L 
is the length of the element. Thus the responsivity of the system is given by 

R = = -L (7) 

This is independent of the loss factors, but the sensitivity is not. The sensitivity is given by 



, 2 1 dPl 1 dP? 

da s = mif = vif (8) 



The fractional power uncertainty from averaging down shotnoise is dP 2 (t) = HujPi/t [15 giving 

^ 2 1 Huj 1 Huj , . 

as ~ L 2 tP[ ~ L2 texp((a G - a L - a s )L)P ^ 



which, for olq = and as <C ol, is optimized with L = 2/aL at 

a\ Huj 
Ttexp(-2)P 



7 9 OL 2 T HUJ , N 

^1=^ 7—7 ONO, ( 10 ) 
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This tells us to make the medium long enough that we see the effect of the absorption, but not so long that the 
shotnoise suffers due to the loss of light. Here we have ignored the effect of spontaneous emission due to any of the 
elements in the cavity. 

We now extend the single pass treatment to an empty cavity, which represents a more realistic practical implemen- 
tation. Advances in mirror coatings and cavity fabrication have reached the state where mirror losses are at the ppm 
level in the NIR. In addition to enabling high finesse cavities, these coatings also have relatively high damage thresh- 
olds which permit large intracavity intensities for measurement sensitivities approaching the photon shotnoise level. 
However, the quoted damage threshold is often significantly above the level where photorefractive processes introduce 
additional noise [16], and 10 kW/cm 2 has been shown to be a reasonable upper limit. Cavity fabrication techniques 
using materials with low thermal coefficient of expansion such as Zerodur [T7] and optimized geometries [18] have also 
rendered high finesse cavities robust to environmental conditions. In addition to mirror constraints, intensity limits 
may also be limited by practical considerations such as power consumption budgets for battery-operated field devices. 

We begin by introducing basic optical cavity terminology following closely the notation of [13] which is used 
throughout. Consider the intensity losses within a Fabry-Perot cavity of length L in one round trip, with each round 
trip taking time 2L/c. The total loss per round trip is given by 



2Lk 

c 



(ii) 



where n gives the total intensity loss rate. The total empty cavity loss is given as 8c = 2Lnc /c = 80 + 81 + 82, consisting 
of the absorptive loss from the mirrors 80 an d the transmissive loss from one mirror given by 81 = ln(l/i?i) ~ 1 — R\ 
(and a similar term for the other mirror with £2), where we have assumed the mirror reflectivity R\ to be close to 
unity. We will assume negligible absorption loss, 8q = and a symmetric cavity 8\ = 82, so that 8c = 2Lkc/ c = 25i. 
The total sample loss is given by 8s — 2Lns/c = 2Las- The total additional absorber loss (or gain) is 8l = 2Lkl/c = 
2LaL, which we allow to be intensity dependent due to saturation. The free spectral range (axial frequency spacing 
between modes) is Af ax = c/2L. The effective finesse of the cavity is F cav = 2ir/8, and the linewidth of the cavity 
(equivalent to the loss rate) is A/ cav = k/2tt = Af ax /F cav . 

We introduce the dynamical equation for the unit less internal cavity field E, where the intensity / = E 2 is equal 
to the number of photons in the cavity mode. We also add a driving term Eq to give [19] 



dE 
~dt 



«(/) 



E + E 



The intensity equation is then 



f f = 2*f = - K (I)I + 2^o 



(12) 



(13) 



In steady state, when all losses are due to the cavity (k = Kc) we require the input power (P ex Eq) to equal the 
output power (HlokcI) which gives Pq = 4HloEq/kc- When the sample and other losses are included we have the 
steady state condition, dl/dt = 0, for the output power 



Pi = fkOKcI 



Po 



4,51 

s(iy 



Po 



(14) 



which can lead to bistability because of the absorption dependence on / through 6(1). 
Ignoring intensity dependence of the absorption, the output power is 



Pi 



461 



(25i + 2La s 



Thus the responsivity is given by 



R = 



6 L f 



Po- 



dPx/das _ -4L _ -2L 
Pi ~~~~5T 



= -L 



2F 



(15) 



(16) 



where we assume 8l = 0, and 28i ^> 2Las- The responsivity is enhanced over the single pass result by a factor 
proportional to the finesse. Note that if the gain is increased past threshold (so 8 < 0), a steady state no longer exists. 
In practice the saturation of the gain limits the laser output power above threshold and is treated below. 
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The uncertainty in as due to shot noise is given by 

2 _ dP[ _ J^dPl_ S 2 hcu 
s ~ (dPi/das) 2 ~W^f~ WL^m 
= (2ft + S L f hw 
16L 2 M£ 2 P 

where we have assumed as is small. This is minimized for Si = Sl/2 = az,L at 



(17) 



da| = = -luT (18) 

Again the optimal output power is not much smaller than the input, with Pi = Pq/4. Note that the sensitivity is 



only a factor of order unity from the optimized single pass case Eq. (10). The difference however is that the cavity 
measurement is much more compact: we can have the same sensitivity for a total amount of trace gas that is l/P C av 
times smaller, in a container that is l/P C av times smaller. 

We next consider the propagation of input 'relative intensity noise ' defined as RIN = dPg/Pg, 

= (A? = (W " ^ WN ' (19) 

We see that the cavity provides robustness to input intensity noise by a factor of the finesse squared in the variance. 
Note that input fluctuations at frequencies larger than the cavity linewidth will be low-pass filtered by the cavity as 
well. 

We also note that if there are technical loss (or gain) fluctuations that appear symmetrically to the absorber, then 
they cannot be distinguished, hence the error would propagate trivially as da 2 s = da\ for loss. Thus, to benefit from 
inserting a medium into the cavity, the loss stability needs to be better than the desired estimation sensitivity. 

In quadrature sum, the form of the uncertainty as a function of time (ignoring non- mirror loss, gain/loss fluctuations, 
and making S\ as small as possible) can be written as 

da 2 = ^!f^-+RIN^) (20) 
s 4c 2 \l + K C t J y J 

where / is the internal photon number and we have used Pi = hwncl and dP 2 = HlukcPi/(1 + Kct)- This form for 
the shotnoise results from filtering the detector output at the bandwidth of the cavity, as fluctuations of the absorber 
cannot be observed at higher frequencies. Note that for CRDS, the RIN component due to input laser intensity 
fluctuations largely disappears and, in the absence of any other noise sources, we are left with the shotnoise limit. 



IV. STOCHASTIC QUANTUM DYNAMICS FOR NONLINEAR METHODS 

So far we have only included optical shotnoise as a fundamental noise source. We have ignored spontaneous 
emission from and saturation of the trace gas, which is appropriate for small amounts of trace gas having broad 
molecular transitions with small cross sections and large saturation intensities. When additional scattering media 
is included however, we must include both saturation and spontaneous emission. To treat these noise sources we 
build on the simplified treatment of Eq. (12) and use a stochastic differential equation (SDE) for the field, with the 
random Wiener increment representing the spontaneous emission (see [20j [21] for comprehensive physical treatments 
of stochastic methods). This treatment for gain within a cavity follows many previous treatments [9j[l2j[T4], however 
much of this section builds upon the analysis of Ref. [12] in particular. 

The SDE with saturable gain is given by 

dEi(t) = Q ( yTW/Is ~ K ') Ei + ^ 0,i ) dt + V / W»(*) (21) 

where the intensity (photon number) is / = \E\ 2 = E 2 + kq = <^g c ? Is is the saturation intensity of the gain 
medium, n' = asc + S\c/L, Eq^ is related to the input power as before, and the Wiener increment dW{ is a Gaussian 
random variable of mean zero and satisfying {dWi{t)dWi{t')) = dtS(t — t'). This SDE is often recognized as a form 
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of the nonlinear van der Pol equation. The SDE's for the fields can be recast into SDE's for the intensity and phase 
using the change of variables 



E 1 = v7cos((/>) 
E 2 = V7sin(» 
dld(j> = 2dE 1 dE 2 , 

where the latter comes from the appropriate Jacobian, and the stochastic change of variables rule [20] to give 



(22) 
(23) 
(24) 



dl(t) = 



d(j>{t) 



Y^JJf- -«')/ + 2Vl(E 0>1 cos(0) + E 0)2 sin($) + 4Q(J) ) dt 



^8Q(I)IdW!(t) 

^ (Eo,i sm(<j>) - E 0<2 cos($) dt + \j^^-dW ck 



(25) 
(26) 



Although generally the spontaneous emission factor Q is a function of /, it becomes independent of / if there is a 
dominant external pumping process, as in a laser. Then the form of Q is related to the gain via Q — Ikq, which can 
be calculated by equating the laser cavity photon number with the atomic number distribution in thermal equilibrium 

Hi. 

These Langevin equations are cast into the form of a Fokker-Planck equation describing the evolution of the 
probability distribution as follows [9 



dEi 



(27) 



where A4 = \ ( i+jffl/i &i + Eo,i and the two noise sources are assumed to be independent. The steady state 

equation is 



dp(E u E 2 ,t) 



dt 



EdJi 



dEi 



2=1 

Ji ~ I Ai 



dQ(I) 
dEi 



p(E u E 2 ,t) + Q(I) 



dp(E 1 ,E 2 ,t) 
dEi 



(28) 



(29) 



where Ji represent the probability currents of the system which must conserve total probability. If Q is independent 
of /, we solve this by setting the probability currents equal to zero, Ji = and integrating by separation of variables 
[9]. In the more general case where Q is a function of /, however, d(A\/Q(I)) / dE 2 7^ d{A 2 /Q{I))/ dE\, and this 
solution technique is inapplicable. (When considering optical bistability, we will let Q{I) = QqI /{I + I s ).) 
Assuming constant Q, Eq^ 2 = 0, and Eq^ = Eq, we have 



Ps(E u E 2 ) 

Pa{I,<t>) 
Ps(I) 



exp 
exp 

r 2 



EqE\ 



(eI + e 2 2 ) + - gIs ^ " 1 ■ ^ 



log(l + {El + E 2 2 )/I s ) \/Q 



EoVlcos^) -"Ll+^L i og( i + I/Ig) ) iq 



/ Ps(I,' 
JO 

exp [(-«'! + k g I s log(l + I /I s )) /4Q] I M bo(EoVi/Q) 



(30) 
(31) 

(32) 



where we have used the definition of the zero-order modified Bessel function 

Imbo(x) = ~ exp cos ((/>)]<# 
n Jo 

2 

7T 



2tt 



exp[a?(l - <^ 2 / 2 )]# 
8 exp(x) 



7T X X 



(33) 
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FIG. 1: (A) Normalized responsivity R as a function of drive field and gain. (B) R versus gain for zero (non-zero) drive field. 
Green solid (red dashed) line refers to the cross-sections plotted in (A) . The dotted black lines diverging near saturation give the 
zero-drive responsivity without the spontaneous emission limit {c/2^'). (C) Intensity as a function of gain for zero drive field 
(green, solid line, sharp kink), and non-zero drive field (red, dashed line, smoothed kink). Parameters are 770 = 10 6 , Si = 10 -5 , 
L = 1 m. 



and x ^> 1 in the approximation. 
Now we calculate the moments 

{n . (34, 

J p a (l,<j>)dl 

and their derivatives with respect to n' which are 

d(I n ) = J °°(-I/4Q)I n p s (I)dI j™I n Ps (I)dI^(I/±Q) Pa {I)dI 

dn> $~p s {i)di {$™p s {i)diy 

= -«/" +1 )-<I")(/))/4Q (35) 
Thus, the responsivity is given by 

_ d(I)/d K ' dK< _ -c (7 2 )-(/) 2 

K - (I) das ~ 4Q (I) ■ [db) 

Note that when the drive is very large, and the system is above saturation, the responsivity returns to that of the 
empty cavity. 

At this point we numerically analyze the entire phase space with the drive and gain, and kg, as the free 



parameters. For the near threshold case, we simulate the responsivity from Eq. (36) using the distribution in Eq. 
( [32] ) . The result is shown in Fig. [l] where the responsivity peaks just above threshold. The threshold curve in Fig. [T]C 
is shown to have a nonlinear kink at this point. But because of non-zero spontaneous emission, the curve is smoothed 
and the responsivity becomes non-infinite, which is the primary result of [12] . The drive term further smooths out the 
curves, reducing the responsivity accordingly, as is shown below. (See Ref. [22] for an exploration of the dependence 
of the threshold point on the drive field.) 

We now consider three different regions of the phase space parameterized by gain and drive: i) laser above threshold 
without drive, ii) laser below threshold with drive, and iii) absorptive bistability with drive. We compare each of 
these schemes to the empty cavity benchmark while considering both quantum and technical limits. 



A. Gain without drive 



Previous demonstrations of intracavity gain medium for absorption spectroscopy have utilized a broadband gain 
medium placed inside the cavity mode. These have included dye lasers, Ti-Sapphire, fiber lasers and diode lasers, with 
the dye lasers demonstrating the best sensitivity [TT]. For compact devices Helium Neon lasers are attractive, specif- 
ically in leveraging existing hardware from a commercial ring laser gyroscope, where the gain medium is integrated 
into a high finesse cavity. As will be discussed below, increased sensitivity is obtained for large laser threshold photon 
number and therefore a tapered amplifier gain medium with saturation intensities at the 1 W level can be considered. 
Using semiconductor lasers may also facilitate gain stabilization because of the phase dependence on carrier number 
[23] . However, careful mode matching of the beam to the micron sized active region of the laser without introducing 
excessive additional loss is required [24 . Nevertheless, advances in high power compact lasers have rendered robust 
intracavity laser absorption measurements encouraging. 

In this section, we keep the gain just sub-threshold such that kq ~ ft', but kq < k! . The gain is also kept low 
enough that the overall intensity is below saturation of the gain element. This section essentially reproduces the 
results of [12 , but with an emphasis on surpassing the empty cavity performance by a given factor in a given amount 
of time in the presence of additional technical noise. 



We consider the limit of / <C i s , where the SDE of Eq. (21) becomes 



dEi = 



KG 



2 \i + \E\yi s 
1 

2< 



k'J Ei + E 0ti j dt + y/2QdW(t) 

(k g (1 - \E\ 2 /I S ) - K ')Ei + £ ,iJ dt + y/2QdW{t) 
= ((k g /2 - k'/2 - \E\ 2 (K G /2I s ))Ei + £ 0) i) dt + y/2QdW{t) 



= ((V -f3'\E\ 2 )Ei + E ^dt+^dW(t) 
and we have translated into the notation of [12] using 



7' = 

a' = 



(«G - «')/2 
kg 
2I S 
q' = Q 

We write the above steady state distribution 

p. (J) = exp [(y//2 - /?'I 2 /4) / q '] I MB0 (E Vl /q') 
-(/ - a) 2 /4\ I M bo(E Ji/^/^) 



Ps(I) 



exp 



where / = iy /3' /q' and a = j f /\Z/3'q'. Without the drive field, the intensity moments are 

- 2exp(-aV4) 
W V7F(l + erf(o/2)) 

~ 2 2aexp(— a 2 /4) 



4exp(-a 2 /4) 
/5r(l + erf(a/2)) ~ tt (1 + erf(a/2)) 2 



(37) 

(38) 
(39) 

(40) 

(41) 
(42) 

(43) 
(44) 



At threshold (a = 0), {I) 
threshold is defined as 



2 



1.128 and (A/ 2 




0.727. The intracavity photon number (intensity) at 



(45) 



for a fully inverted medium with Q = Ikq- The below saturation approximation (I) <C i s , implies a <C 770 and 
(J) <C (A/ 2 ). Furthermore, if we assume the below saturation limit, the above threshold limit a ^> 1, and 770 ^> 1, we 
have 



KG 



4a 



(46) 
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Thus the following limits are equivalent 



K 

7t 



a 




(48) 


y 


< K 


(49) 


(i) 


<C Is 


(50) 



In the limit of large a ^> 1, the normalized quantities are (J) = a and (A/ ) = 2, which give 



(J) = £1770^=27,7'/*' (51) 
(A/ 2 ) = 2 % 2 J=8J S . (52) 



The responsivity calculated from these moments is essentially the same as in Eq. ( 16 ) 

d(I)/das —c 



R (I) 2V 



remembering that 7' depends on ft' which depends on as- The enhancement factor over the empty cavity, Eq. (16), 
is Si/ 5 which can be made infinite at a = 0. However, this result is invalid unless a ^> 1. When the full equation is 
used, the responsivity, as displayed in Fig. [Tj3, peaks at 

R = d -^ = -^(°^vi:) (54) 



(I) Si 

for a = 1.28 (the normalized function d ^^ da i s maximized at 0.346 for a = 1.28) [12]. Thus, the responsivity 

enhancement factor over an equivalent cavity without a gain element is nearly vTs • 

We now determine the sensitivity da 2 s . We perform a linearization of the dynamics which allows us to easily apply 
standard estimation theory techniques. We start with the full nonlinear problem with Q constant as given in Eq. 
(21). For the linearization, we assume E\ = Ei + fri, with \bi\ <C E%, where Ei solves the equation without noise, 

°'l(TTwfElW,- K ') E ' + ^ 1551 

We also ignore all b\ terms and assume £^0,2 = 0, E 2 = and consider only b\. These approximations give 
dE x = d{E x +61) = dh 

= 4((irlk)((Tflfe)- 1 ) + ' t ') Mt+ ^ <iW ' (<) (56) 

In the far below saturation limit E\jl s <C 1 and the dynamics are 

1 



dh(t) = ^(K G (l-ZEl/I s )-K , )b 1 dt+^dW i (t) (57) 

With no driving field, E\ « I s (k>g — k')/kg (need above threshold kq > k') and 

dbi(t) = - (k g - K f ) bxdt + y/2QdWi(i) 

= -2 7 'Mt + y/2QdWi(t) (58) 

which is an Ornstein-Uhlenbeck process with mean zero and variance (Abf) = 2Q/4 / y / . Because we have assumed 
that the system is above threshold, kq > k' this equation is always damped. Whereas in the full equation the 7' 
term acts as a gain, in the linearized form it turns into a damping. Because I = E\ + 2E\b\ + b\ « E\ + 2E\b\ with 
Ei a constant, this gives 

(A/ 2 )o = 4£ 2 <A6 2 > = 4^ = ^ = 81. (59) 
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using kg = k'. To check the validity of the assumption we also need 



1 « = 2a2 ' (60) 

thus, the linearization is valid in the regime 1 « a « 

The relative intensity uncertainty as a function of time is 

(A/ 2 )ft) = (A/ 2 )o/(/) 2 = M s /(2I s j/n') 2 = 2n' 2 /I s j 2 

(i) 2 i + 2yt/2 i + yt i + yt [ } 

where j' is the effective bandwidth of the fluctuations. Comparing this to the shot noise limit of a coherent state 
within the cavity, which for output Pi = hwK 1 '(I) gives 

(AI 2 )(t) = K^Prin = W)_ = ( , 

(I) 2 1 + l + K't l + K't [ } 

where n' is the cavity linewidth. Below saturation 7' <C ft', the shotnoise variance is less than the fluctuation variance 
above and negligible. This shows how the laser case limits to the empty cavity case as 7' approaches k' . 



We now systematically compare the performance of the total gain estimation uncertainty (from Eq. (53) and Eq. 



(61) with technical noise added in quadrature) to the empty cavity uncertainty of Eq. (20), which are given by 



o, 4 7 /2 /2ft' 2 1 x 
(da\) G = 4- — — + v T (63) 



2 \i s y 2 i + yt 

= ^(srn5« + *) (64) 

where the subscripts G and E denote gain and empty cavity respectively. In the empty cavity case, we have replaced 
the RIN with vt for the sake of generality. 

When comparing the two cases we make the following assumptions. First, we assume cavity dimensions are the 
same but set n' G ^> n' E by an order of magnitude to account for extra loss in the cavity with gain due to mode- 
matching through the gain element. We assume that there are no technical gain or loss fluctuations that would affect 
the nonlinear case more severely, and we also assume that vt (due to an unspecified technical noise source) is the 
same for both. We optimize the internal intensity of the gain case between the ranges of validity of the expressions 
used vTs <C (/) <C i s , but also limit it by the same mirror threshold. The mirror threshold intensity Im in units of 
photon number (and Im,u m regular intensity units) is Im = Im,u]^; an d we restrict the gain case by (I) < Im- We 
also set the empty cavity intensity at this limit, Ie = Im- 

For short times, fluctuations in the gain case are worse than shotnoise, but with vt = and an equivalent number 
of internal photons, the long time (averaging down as 1/t) limit is approximately the same for the two cases. 

Fig. [2] shows the estimation performance as a function of time. For the gain case, the uncertainty starts at the 
variance level of the cavity dynamics, integrates down at times greater than the inverse of the fluctuation rate until 
it hits the technical noise floor set by vt- (The shape of the empty cavity case is the same, but relatively flattened 
because of the high value of the internal intensity, as vt ^> 1/Ie-) The maximum improvement occurs at long times 

We find that one should make 7' small but this slows the fluctuation rate and reaching this limit takes too long. 



Instead, we consider the case where we pick a specific time (one second) and optimize Eq. (63) at that time over the 
variable 7' . If the optimal 7' falls outside the allowed range n' G / \/% « 7' « k! g (or the intensity goes greater than 
mirror threshold), we fix it at the closest bound. 

Using this approach, we plot (dot's) c(t) m Fig. [3^ which is optimized at 7' corresponding to the internal intensity 
shown in Fig. For sufficiently large vt, the gain case can significantly outperform the empty cavity limit of 
(do? s ) E(t) , although both clearly get worse in an absolute sense as vt increases. In this sense, the gain technique is 
more robust to the technical noise than the empty cavity approach. 

We analyze the robustness improvement factor as follows. At some time t x we desire the sensitivity improvement 
with gain to be some factor x better than the no gain performance. We define the time that the gain- sensitivity levels 
out as to — 2k'g/ I s Y 3v T, while the time the empty cavity sensitivity levels out is ts — 2/(IeVt^)- The intersection 
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FIG. 2: Estimation performance for a cavity with gain (green solid line) and without gain (red dashed line). The gain case 
eventually integrates down past the no-gain case in the presence of technical fluctuations due to the enhanced responsivity. 
The gain case is optimized over a for peak performance at one second (vertical dashed line), and two curves with a above 
and below the optimum value are also shown (green dashed-dotted line). Parameters are 770 = 10 6 , S liE = 10" 5 , S liG = 10" 4 , 
A 1.064 /xm, A tt(1 mm) 2 , L = 1 m, I M ,u 10 kW/cm 2 , v T 10" 9 . 



time for when sensitivity with gain becomes better than that with an empty cavity is given by setting the empty case 
long time limit to the gain case intermediate time and is 



c 2 i s yt c 4c 2 

7f2 



vt (66) 



tc = ,2V"? • ( 67 ) 



We get the proper constraint by setting t x = xt c < t G - Note that this does not depend on Ie because we assume 
the empty cavity performance hits the limit due to vt earlier than the crossover time, i.e. ts < t c . We can rearrange 

V = 25 ^ X (68) 

In other words, for a fixed vt, for t x to be reasonably small and \ reasonably large, we need 7' to be reasonably large 
and away from threshold (but still less than k/ g ). We also need x^c < to or 

tn 16^1 

X< Tc =^- (69) 

Numerically, this constraint describes the performance in the large vt limit. In Fig. [3j the two sensitivities cross each 
other (x = 1) when the technical noise variance is equal to 

Below the critical value of vt, the empty cavity wins, and above the critical value the nonlinear method is preferred 
with an optimal photon number that decreases as vt increases. Here the sensitivity improves by the factor in the 
inequality above until the system approaches threshold (7' ~ n' G /^/T s or a ~ 1) at which point the analysis becomes 
invalid (and does not diverge when made valid). We also note that for small enough vt the empty cavity case levels 
off due to shotnoise. 



B. Gain with drive 



In this section, we demonstrate how spontaneous emission prevents gain added to a cavity with a drive field from 
producing an effectively 'infinite finesse' cavity. We have previously seen in Fig. [TJthat the responsivity with drive 
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FIG. 3: (A) Estimation performance in one second for the empty cavity at the mirror limit intensity (red dashed line), and the 
cavity with gain at optimized a (green solid line) . (B) The green solid line gives the internal photon number corresponding 
to the optimal a, with a bounded by 5 <C a <C r]o {y/Tl <C (I) <C I s , green dashed-dotted lines). The red dashed line gives the 
internal photon number used for the empty cavity case (at the mirror limit, Im,u)- Parameters are the same as in Fig. [2] 



does not do better than the zero-drive case, but it can outperform the empty cavity case in the presence of noise. Here 
we aim to replicate the analysis of the last section with the goal of maximizing the sensitivity over both drive and 
gain for a fixed amount of measurement time. In this case one can fill the cavity with the drive field to optimize the 
internal intensity, instead of relying on gain alone. This may have practical consequences when considering technical 
constraints. 

We follow the same linearization treatment as before, starting with Eq. (21) and Q = Ikq- However, we now 
set 7' = (k' — Ko)j2 > (so loss is slightly greater than gain) and also assume 7' <C n' ~ kq. To ignore the 
saturation term entirely, we require the condition kq! /I s <^ d — — 27'. We can then linearize with E\ = E\ + 61, 
I = Ef = Eg/y 2 , where E 2 Q = E^ and Eg >2 = 0, giving 



dbi(t) = -^ l b 1 dt+^2QdW i {t) (71) 
The variance is (A&f) = 2Q/2Y = 2kg/j'- For the linearized treatment, we need I ^> (A6f), or 

El > 2k g7 ' (72) 

and to ignore the saturation term we need 

k g I/I s < 2i (73) 

E\ < 2I s1 ' 3 /k g (74) 

Note that this is more stringent than / = .Eg/7' 2 < I s because 7'/ K G •C 1. Putting these conditions together we have 

2k g1 ' < Eq < 2I s1 ' 3 /k g (75) 

i > k g /^/T s (76) 

Thus for 7' we have 

kg/V^T < 7' < kg (77) 
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which gives 

y/T a <£l= Elh 12 « I s (78) 
Using I = Eq/j' 2 and 7' = (kc + cas — «g)/2, the responsivity is 

= dJI ^ = {di/das) d -^- = ^ (79) 

I i7 

which since 7' <C ^g, can again be better than the empty cavity by a large factor. The variance is given by 
(A/ 2 ) = AI{b\) = 4(Ei/y 2 )2K G / 7 ' = 8 KG E 2 / 7 ' 3 . 
We now minimize the sensitivity 



v'2 
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(A/ 2 ) 


1 


I 2 1 


f 7'V 2 




1 


El 1 


+ 7'V 2 


4k' 2 


1 


LY 2 1 4 


-yt/2 



Vt 



v t (80) 



where in the last step we maximize Eq at its upper bound and use kq ~ d . This is essentially the same as Eq. (|63b 
with the same limits on 7'. Thus, when optimizing over 7' at a specific time, we observe the same trends as in Fig. [3] 
for large vt there is improvement optimized with 7' small, and as vt decreases, the optimal 7' increases and it gets 
more difficult to beat the empty cavity performance. 



C. Saturable Loss 



In the gain-enhanced example above, the nonlinear kink near laser threshold as represented in Fig. [T] provides 
enhanced performance in the presence of technical noise for a steady-state measurement. Within the same phase 
space, but operating with saturable loss rather than gain, increased responsivity is expected due to optical bistability, 
which occurs when the intracavity intensity is near saturation of an internal absorber, with the below-saturation 
absorber loss being much greater than the empty cavity loss p~3j [14] . 

This system is analogous to a particle within a drive-dependent potential. In the bistable region there exists a 
double-well potential, but if the cavity is driven too high or low, the potential is tilted such that only one local 
minimum exists. In the quantum case, spontaneous emission induced tunneling between the wells occurs and the 
system equilibrates over long timescales with a single- valued mean for all drive fields. The slope of the output versus 
input curve becomes large in the middle of the bistable region where the system spends an equal amount of time 
in each potential well. Thus the variance is also large at this point corresponding to a large peak in responsivity 
[T4] . Note, however, that the tunneling rate becomes slow near this point, thus there is again a trade-off between 
fluctuation timescales and responsivity. 

Implementing absorption spectroscopy with saturable loss requires temperature stabilizing a cavity with an internal 
saturable absorber such that the addition of the trace absorber kicks the intensity from one stable point to another in 
the bistable hysteresis curve. In this section, we show, as in previous examples, that spontaneous emission limits the 
fundamental sensitivity from surpassing that achieved with CRDS, but once again, this method can make the system 
robust and the effective sensitivity superior to CRDS in the presence of technical noise. 

We begin with the classical optical bistability equations and make the absorption intensity dependent through the 
relationship 

«(0 - YTJJT,' < 81 > 



From Eq. ( 14 ) we have 
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FIG. 4: The 'sweep- up' modulation technique simulated for a 10 -8 cm -1 differential absorption signal and kl — 12^c (greater 
than the 8k,c needed for bistability ) . (A) The output power (for the two cavities, red and blue) is shown versus the common 
input power, which is swept over the steady-state bistability curves (dashed lines). The top gray line represents input equal to 
output. (B) The two outputs (red and blue solid lines) as a function of time for the given input ramp (black dashed-dotted 
line). The differential output power (magenta dashed line) reveals the degree of differential absorption via the switching times. 
Parameters are the same as in Fig. [2] but with an absorber having saturation intensity I s , u — 1 W/cm 2 . 



Plotting this equation gives the s-shaped optical bistability (hysteresis) curve. The two internal intensities at which 
the slope is infinite are given by setting the derivative dPo jdl = and solving which gives 



T kl~ 2(kc + ca s ) ± >/kl(kl ~ 8(^c + ca s )) 
2(k c + ca s ) 



(83) 



We define a unitless parameter x <C 1 such that kl = 8(k,c + ca s)(l + x 2 ), and I± ~ I s {3 ± Ax). We see that the 
bistability occurs for kl > 8^c at values greater than three times the saturation intensity. The resulting bistability 
curve exhibiting hysteresis is shown in Fig. [4] for different levels of trace gas absorption. 

Instead of repeating the steady-state treatment as in the examples above, we now consider a modulation based 
approach for several reasons. First, the steady state Fokker-Planck solution cannot be integrated as above because 
of the intensity dependence of the spontaneous emission parameter Q(I) = QqI/(I + I s ). Second, it is natural to 
compare to the CRDS standard which is also a modulation approach that is robust to several technical noise sources. 
Third, we stress analogies to other modulated nonlinear metrology systems, such as fluxgate magnetometers, where 
the signal is sensed via a shift in a measured switching time. 

We consider two potential modulation schemes. First, we analyze a 'sweep- up' scheme where the drive is ramped up 
and the difference in switching times measured. Second, we compare ring-down schemes with and without saturable 
absorbers, and show that including the saturable absorbers allows the system to outperform empty cavity CRDS for 



certain technical noise models. In both cases, we simulate single stochastic trajectories using the SDE of Eq. (21), 
but with = -k l < and Q(I) = QqI/(I + I s ). 
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FIG. 5: As the differential absorption level is decreased the signal becomes obscured by spontaneous emission induced timing 
jitter. For a signal of 10 -10 cm -1 (magenta) the differential blip is dominantly in the same direction (negative), yet when the 
signal is decreased to 10 -12 cm -1 (black) the order of hopping, hence the sign of the differential, becomes random, thus giving 
the noisefloor. Parameters are the same as in Fig. [4] 

1. Sweep-up with saturable absorbers 

We define the sweep-up modulation scheme as follows: we consider a measurement scheme with two cavities mounted 
in parallel on a monolithic block, and admit the trace absorber into one cavity but not the other. Alternatively, with 
a broadband saturable absorber but a narrow sample absorber one can use two frequency beams, one resonant and 
one far detuned from the sample. Both cavities or frequencies are driven with the same laser source dithered across 
the entire dual- valued part of the S-shaped response curve shown in Fig. [4] As the drive field is swept to increasing 
values, the intracavity intensity will jump to the higher branch at different times due to the absorption differential 
from the trace gas. 

In this modulation scheme, the sweep rate is constrained to be faster than the tunneling time between stable points 
such that the intensity can be deterministically swept to the jump point without the internal intensity randomly 
switching states. But the sweep rate is slow compared to the ring-down rate of the empty cavity kc, or else the 
jump will be smooth. Simulation results for this scheme are shown in Fig. |4j where the input intensity is swept 
over the bist ability region at a moderate speed. For the large differential absorption used in the plot, the difference 
in switching times is represented clearly in the difference of output intensities. The differential intensity for smaller 
signals is shown in Fig. [5j If the differential absorption is large enough, the differential signal always has the same 
sign. But as the absorption signal falls below a certain value, the order of the switching times become randomized due 
to spontaneous emission induced tunneling. This scheme is analogous to the fluxgate magnetometer design, wherein 
a small field is sensed via the differential switching times of two oppositely driven solenoids containing ferromagnets 

2. Ring-down with saturable absorbers 

Another modulation technique involves loading the cavities with an intensity much greater than the saturation 
intensity and monitoring the differential intensity decay or ring-down curve as in CRDS. As opposed to CRDS, the 
evolution with saturable absorption displays two different timescales corresponding to the loss rates both above and 
below saturation. The system starts decaying at a rate kc as if the absorber were absent, but once the intensity hits 
the saturation level, the absorber turns on and the decay becomes much larger, kc + ftc? as shown in Fig. [6^L The 
transition time between these two regimes depends on the internal trace absorber, and this time can be estimated by 
recording when the intensity drops below a certain sub-saturation threshold as compared to the other cavity. 

Again the timing performance will be limited by the spontaneous emission induced random tunneling. We determine 
the absorption estimation error from the timing error as follows. The intensity starts at imit? an d decreases past I s 
at time t from which we can measure the absorption from as = — (log(Ii n i t / 1 s )/ct) — k>c/c. When we are ramping up 
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FIG. 6: (A) The dashed red curves show empty cavity ringdown as a function of time (in semilog scale) with 'measurement 
noise' separation of g — 0.5^c using the model described in the text. The same measurement separation is shown for the 
saturable CRDS in the solid green curves. The absorption signal is the same in each case and, in the saturable case, estimated 
from the times that the curves pass through a threshold intensity (green dashed-dotted curve at is/10). (B) The performance 
from fitting regular CRDS (red dashed) and timing estimation with saturable CRDS (green solid). Because the green curve 
timing is relatively immune to the multiplicative measurement noise, it performs better for large amounts of this noise where 
the CRDS estimation is limited by the indistinguishable measurement noise. For low amounts of noise, the red curve hits 
the shotnoise limit and the green curve hits the limit from Eq. (84). These approach each other as the saturation intensity 
approaches the mirror limit intensity. We use the same cavity parameters as in Fig. [H] but with kl — lOO^c, Is,u — 1 W/cm 2 , 
and initial intensity at the mirror limit of Ii n it,u = 10 kW/cm 2 . 



or down near I s at rate Kc the variation dl will trigger the jump process in a time window of width dt = dIjI s Kc- 
Thus the total sensitivity on a single shot will be roughly 



log ( ^ ) ^ « log ( ^ ) — ^ « — (84) 



dt . ( linit A K C dl K C 1 



I s J ct 2 y I s J C I s 

where we have used t ~ 1/kc an d dl « \/T Sl log(/i n it// s ) ~ 1- This is roughly the same as the fundamental limit of 



Eq. (64), but with I s instead of the Because I s cannot surpass Ie this limit is generally worse, but we see that 
saturable absorbers with high I s are preferable. This same argument can be applied to the sweep- up case of the last 
section as well. 

We next compare the robustness of this saturable ring-down method to empty cavity CRDS for a particular noise 
model. We propose a certain multiplicative measurement noise model that takes the output light and multiplies it by 
a factor exp(g't), where 'gain drift' g' in inverse time units has mean zero but finite variance and is constant within a 
single shot (i.e., (g' 2 ) ^ 0, (g') = 0). This model represents noise due to a time-dependent gain in the measurement 
process which could, for example, be due to actual gain fluctuations in the detection or alignment instabilities. This 
is clearly seen when the exponential is expanded to give 1 + g't. However, the exponential form of the model was 
chosen such that this noise process is fundamentally indistinguishable from the cavity decay in CRDS. In other words, 
the noise takes kq — > kc — g' and thus no fitting procedure can distinguish kq from g' in a single-shot. Thus, the 
estimation error is equal to (g /2 ) when this is greater than the shotnoise limit. In contrast, the nonlinear ring-down 
method displays a kink which is not shifted substantially in time when this noise model is applied. The threshold 
value from which the timing is taken should be below the switching point where the intensity becomes less than 
saturation I < I s . Past this point the loss curve is sharper vertically by a factor of kl/kc ^> 1, hence the timing is 
robust to such time-dependent gain errors by this same factor. The resulting robustness of the nonlinear estimator is 
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compared to the non-robust CRDS in Fig. [6j The timing estimator is robust enough to outperform CRDS for large 
values of (g f2 ) but the error does begin to rise at extremely large values of this noise. 

This multiplicative measurement noise model is generally representative of any time-dependent measurement gain 
during the time of the ring-down. One could propose another noise model that moves the kink location in the nonlinear 
case as well, but such a model is practically less likely for slow noise because it relies on having at least a quadratic 
noise signal in time, requiring a higher order expansion of the noise source. We stress that this improved robustness is 
only for a particular technical noise model on the measurement side. The CRDS method remains relatively immune 
to initial intensity fluctuations, whereas the saturable-absorber method does not. 

V. CONCLUSIONS 

In this paper, we have analyzed nonlinear methods for high-sensitivity, intracavity absorption spectroscopy. Al- 
though inherent quantum noise mechanisms, such as spontaneous emission and shot noise, prevent the fundamental 
limits of nonlinear approaches from surpassing those of linear methods (e.g., empty cavity ringdown spectroscopy), 
we have shown that they present an advantageous robustness to certain kinds of technical noise. This enhancement is 
due to nonlinearities of the system response which can be a function of many experimentally convenient parameters, 
such as drive, gain, or time. 

Implementing these nonlinear techniques involves additional experimental costs and challenges, which may ulti- 
mately introduce larger technical noise sources. The gain and saturable-absorber systems will have some degree of 
technical fluctuations in the medium (e.g., due to thermal atom number fluctuations) and also will likely reduce 
the cavity finesse due to insertion loss. However, specific implementations may warrant the effort when considering 
commercial, field- deploy able devices where environmental noise sources may be significant. 
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